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ABSTRACT 


This report summarizes the testing of the various magnetic 
field models against the available World Magnetic Survey data 
and describes the method by which the first International 
Geomagnetic Reference Field [XGRF (10/68) was derived. The 
IGRF (10/68) was composed of contributions of the field models 
derived by: Goddard Space Flight Center, Air Force Camoridge 

Research Laboratories, Royal Greenwich Observatory, Institute of 
Terrestrial Magnetism and Radiowave Propagation (IZMIRAN), and 
the U, S. Coast and Geodetic Survey. 

IGRF (10/68) is a set of 80 internal spherical harmonic 
coefficients and their first time derivatives epoch 1965.0 
referenced to a sphere 6371.2 km radius. The root -mean-square 
residuals to surface and airborne magnetic survey data taken 
during the interval 1961 through 1965 average about 200y. The 
rms deviations from selected COSMOS-49 (1964.7) and POGO (1965.8 
1967.9) total field satellite observations range from 30 to 60y. 


Introduction 


This is to summarize some of the computations that took place 
at the IAGA Symposium in Washington, D. C., October 22-25, 1968 
that led to the resolution by the Working Group on the Analysis of 
the Geomagnetic Field (Reporter, A, J. Zmuda) to propose a particular 
International Geomagnetic Reference Field. The basic requirements 
established by Dr. Zmuda following the discussion at previous 
meetings was that the IGRF would consist of no more than 80 spherical 
harmonic coefficients of internal origin epoch 1965.0, with each 
having a first time derivative. These coefficients were to be true 
spherical harmonics in describing the field as opposed to those 
M quasi-spherical" coefficients resulting from derivations neglecting 
the oblateness of the earth. Further, only sets of coefficients 
submitted to the Working Group on or prior to March 15, 1968 were to 
be considered. 

These sets of spherical harmonic coefficients are given in Table 1. 
They are each updated to 1965.0, and limited to an n* (maximum degree 
n and order m) of eight. Of the sets given, all except (g) and (h) 
meet the requirement of taking into account the oblateness of the earth 
in their derivation. Most of the field descriptions are to appear in 
the WMS volume to be published. However, a few have been published 
separately as follows: 
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Table 

Field Model 

Reference 

*r 

la 

GSFC(l2/66) 

Cain et al. * 1967 


4 

USC&GS 

Hurwitz et al., 1966 

9 

Ih 

RGO-l(LME) 

Leaton et al., 1965 


Test Data 

Although no explicit formula was agreed upon prior to the 
meeting for the derivation of an IGRF, there was an understanding 
that the considerations would need to be somehow based on the 
correspondence between the field components predicted by the 
proposed models and the available survey data. 

Since the epoch of this IGRF was to be 196^, an arbitrary data 
cutoff time of 1961 was chosen so the results would not be too 
heavily weighted by observations prior to 1965» Testing was done 
on all data available since that date. These were divided into the 
major categories as follows: 

a) surface magnetic observatory annual means 1961-1967 

b) surface magnetic surveys. This category includes land 
surveys, repeat stations, shipboard and ship-towed 
observations . 

c) aeromagnetic survey of Japan (1965) f Nagata , 1966] 

d) aeromagnetic survey of Canada (1961-1963) 

e) aeromagnetic survey of Scandinavia (1965) 

f) project MAGNET worldwide (principally oceanic) airborne 

survey (1961-1966) f USNOO , 1965] 
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g) OGO-2 data as available during magnetically quiet 
intervals October, 1965 - September, 1967 

h) 060-4 data during magnetically quiet Intervals 
July - December, 1967 

i) 1964-83c observations - 1964-1965 f Zmuda et al, , 1968] 

j) COSMOS -49 observations* - 1964.8 

k) Other airborne (towed proton magnetometer data) 

All of the non-satellite data were obtained from the file 
prepared by the Geomagnetic Division of the U,S. Coast and Geodetic 
Survey (E.Fabiano and S. Gain, WMS Volume, 1969), This file contained the 
contributions from many separate organizations a^d survey groups 
and is constantly updated as new observations are submitted, This 
file was edited by rejecting those observations deviating from the 
GSFC (12/66) model (using n* = 10) by more than lOOOy. This 
procedure was used to eliminate the highly anomalous data beyond 
about five times the rms (root -mean-square) deviation. Since all 
models were truncated to n* « 8 for testing, it gave no particular 
advantage to GSFC(12/66). This model was sin-re it fitted the 

data set best, hence eliminating the smallest fraction of data. 

The amount rejected is seen to be small as given in Table 2: 
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TABLE 2 

NON-SATELLITE DATA ELIMINATED FOR AC > lOOOy 


Component 



Data Type 

Observations* 

Data 

Reiected 




* No 

% 

a) 

Observatory 

1984 

34 

1.7 

b) 

Surface 

22425 

204 

.9 

o) 

Japanese Air 

1461 

6 

.4 

d> 

Canadian Air 

9470 

27 

.3 

• ) 

Scandinavian Air 

6973 

1 

.01 

f) 

Project MAGNET 

104228 

401 

. 4 

k) 

Other Air 

1763 

9 

.5 

In this and ensuing 

discussion a value 

of D, I, H, Z, or F 


is counted as one observation even though other values may have 
been measured at the same time and location. 

The OGO-2 and 0G0-4 data (sampled every 30 seconds or at a 
spacing of approximately 200 km) were initially selected from 
periods of time for which Kp = 0. They were then fit with a special 
model listed in Table 3 [POGO(10/68) ] employing 143 internal coeffi- 
cients and their first time derivatives. The distribution of 
deviations of the data from this fit was as follows: 

|af|y 0 10 20 30 40 50 60 70 100 200 600 Total 

Obs. 27646 4218 589 141 23 26 6 2 9 4 32664 
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Since the distribution indicated that the 15 observations 
over 70y were Hkely anomalous, they were rejected and the resulting 
rms deviation computed to be 7y, The remaining 32649 observations 
were included in the testing. 

The COSMOS-49 data were similarly treated by fitting with a 
special function and eliminating those data that deviated signifi- 
cantly from the rest. The data were prepared by the U. S. Coast and 
Geodetic Survey from the catalog published by IZMIRAN ( Dolginov et al . . 
1967). These were sorted into time order and each fourth observation 
fit with a series or 99 spherical harmonic coefficients. Data 
exceeding lOOy from the fitting surface were rejected in the coefficient 
determination. The distribution of residuals from this model, labelled 
COSMOS (9/68) , is as follows: 

| AF | y 0 10 20 30 40 50 60 70 80 90 100 Total 

Obs. 1853 1243 648 271 93 41 '23 18 19 15 138 4362 

The use of every fourth observation in the fit is adequate for 
these purposes since each orbit then contains about 10 observations 
for the shortest wavelength used of the fitting function (n* = 9 
corresponds to 360/9 = 40°). Since the RMS deviation of these data 
from the COSMOS (9/68) field was 21y, the selection used for model 
testing were those deviating less than 60y, a total of 16554 from 
the approximately 18,000 originally available. 

The 1964-83c observations entered the testing unedited except 
for the rejection of one spurious point that gave a |AF| > lOOOy. 
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Test Results 

The various models were tested against the data sets both with 
the limitation oi : 80 coefficients and also using all coefficients 
if more were available. Table 3 illustrates for the GSFC(1.2/66) 
model the distribution of residuals using the first 80 coefficients 
as well as the full number. Since the surface data were edited with 
this model using a 1000y criteria, there can be no residuals above 
this figure with 120 coefficients. The effect of the truncation is 
seen to increase the rms residuals by 10-20y regardless of their 
magnitude. Using 80 terms has only a small percentage effect on 
the surface data since magnetic anomalies account for a great deal 
of the scatter. The consequence on the satellite data is more 
obvious as seen in the 0G0-2 results. Here the effect is to increase 
the observations in the 50~100y range from 5 to 10% of the total data, 
and to push the number over lOOy from 1 to 3%. 

These distributions were also calculated for each of the other 
test models and the rms values compiled in Table 4. Here the relative 
match of each data set to each model can be readily compared. 

Although for each model there is an improvement with an increase 
in coefficients, the differences are generally smaller for those with 
higher average residuals. 
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Root -Mean-Square Deviations of Test Data From Various Models 
Using n v ' c = 8 and n* = Maximum Degree and Order of Expansions 


% 


oo 

vo 


o o 

H rH 


i 

SI 


Ph 

c/3 

O 


* 

CS 


VO 

ON 


rH 

r*. 

CO 

r*» 

rH 

O 

CM 

ON 

O 

CM 

CM 

VO 

o 

in 

CO 

CO 

CM 

CM 

CM 

H 

CM 





CM rH 

r** r*. 

Cs| CM 


00 00 

in 

CM CM 


On <t 

n 

CM CM 


<t r^. 

co co 

CM CM 


n co 

in in 

CM CM 


o n 

co CM 
co co 


r**. on 

Mf rH 


CO rH 
CO CO 


<r 

On ON 


CO'H 

CM 

CO 

n 

co 

CM 

VO 

ON 

00 

rH 

O 

in 

n 

CM 

VO 

rH 

On 

in 

CO 

CM 

CM 

CM 

CM 

rH 

CM 




vO 

VO 


00 

VO 

oo 

■O CM 

CO rH 

00 H 

00 

ON <J- 

r>* rH 

r*. in 

n* 

■***. 

CO 

O ON 
CM rH 

O ON 
CM rH 

*0 CM 
CM CM 

CO CM 
CM CM 

i**- vo 

rH rH 

O O 
CM CM 

VO VO 


n 

r** 

vO 

00 rH 

<t ON 

CO VO 

O CM 

CO O 

vo r*. 

O I s *. 

oo r*. 

CM 

's 

rH 

rH 

O O 
CM CM 

rH O 
CM CM 

CM rH 
CM CM 

rH 
CM CM 

vo n 

rH rH 

rH rH 
CM CM 

00 I s * 

VO VO 

uo 


ON 


o 

co 

00 

vO 

rH 

CO VO 

<f in 

r*>. on 

ON 00 

CM 

rH 

o 


rH O 

O 00 

CO rH 

CM O 

in co 

CO 

CM 

PH 

CO 

CM CM 

CM rH 

CM CM 

CM CM 

rH rH 

CM 

CM 


ON O 
<fr CO 


*0" r— I 

CO CO 


00 

CM 


X 

00 CO 

6 


_ X 

00 to 

a 


X 

oo rt 
B 


X 

00 to 

B 


» a 

6 


00 to 

B 


x 

00 CO 

B 


X 

00 to 

B 


X 

oo co 
B 


in 


C/D 









O 00 

<t CM 

in vo 

rH CO 

VO CO 

ON VO 

1"* O 

-0 I s * 

VO ON 

VO CO 

% 

<f CO 

CO CM 

!**■ CO 

*0 CO 

ON ON 

co 

<t CO 

ON ON 

rH O 


CJ 

C/3 

P 

CM CM 

CO co 

CM CM 

CM CM 

rH rH 

CM CM 

rH rH 


rH rH 

rH H 


*o 


ON 

00 


CO 

O 

ON 

O 

in 

•o 

ON 

in 

00 

vo 

00 

CM 

ON 


CO 

00 

CO 


oo 

ON 

CM 


CM 

CM 

CM 

CM 

rH 

CM 

rH 



rH 



ON o 
00 00 


«n cm 

oo oo 


o 

00 

VO 

CO CO 

r*. r** 

ON O 

vo in 

ON O 

in cm 

rH 

CO ON 

o r**. 

O ON 

o 


O ON 

o oo 

CO CM 

CM O 

in <t- 

rH O 

m cm 

CO CM 

to 


o 

Pi 

o 

rH 

CM rH 

CM rH 

CM CM 

CM CM 

rH rH 

CM CM 






on in 

CO rH 


VO 

vo 

00 I s *. 

CM O 

VO rH 

in CO 

cm in 

O VO 

00 I s *. 

CM 00 

ON 00 

r-4 



ON 00 

O 00 

CM rH 

CM O 

vO 

O 00 

<t CM 

CO CM 

CO CM 

■VO 


CM 

rH 

rH rH 

CM rH 

CM CM 

CM IN 

rH rH 

CM rH 







X 

00 CO 

B 


< 

O 

rH 

in 

CO 

CM 

1**. 

<t 

o 

CO 

vO 

H 

n 

CM 

n 

Mi- 

f*^ 

CM 

in 

CO 

ON 

in 

< 

ON 

CM 


nt 

ON 

00 

n 

to 

<r 

rH 

O 

rH 

CM 

rH 

ON 

VO 

CO 

VO 

rH 

ON 

CO 



CM 




o 

rH 


rH 

rH 







vH 











H 










a 

g 






>N 




to 

O 






u 




T? 







o 






ON 






CO 

<u 

a 

to 



O 





<D 

to 

to 

c 

U 

t 

CO 




> 

C/D 

a> 

•H 

•H 

o 

C/3 

00 




u 


a 

P 


<y 

Q 

1 

CM 

*o 


<u 

'O 

to 

co 

a 

•n 

X 

<t 

1 

1 


<0 

a 

a- 

c 

CO 

o 

CO 

vo 

o 

o 


fO 

CO 

CO 

CO 

O 

M 

o 

ON 

o 

p 


o 

p 

p 

o 

CO 

Ph 

p 

rH 

o 

o 


00 


CM 


ON 


00 


ON 


% 


Max . value of 



We ighting of XGRF 


It was decided that a weighted average of coefficients would 
provide the best compromise tc an IGRF < With the restriction that 
models to be included should be based on truly spherical coefficients, 
the RGO(LME) and USGGS models were eliminated for inclusion in the 
main field averaging. Since the surface data residuals were so 
greatly influenced by crustal anomalies, it was decided to base the 
weights on the residuals to the satellite data. 

Several different weighting schemes were tried. Generally, 
the exact weights used did not alter the overall results appreciably 
as long as those models fitting the satellite data best had an 
advantage. The POGO(3/68) and AFCRL(ll/67) models were eliminated 
from the considerations since there was another model submitted by 
the same organization. 

After several semi -qualitative arguments and considerations 
that the IGRF should be most useful near 1965.0, the following table 
of relatively weights were agreed upon to be applied as inverse square 
factors in combining the main field terms: 


Model 

a 

GSFC(12/66) 

40 

AFCRL (3/ 68) 

70 

RGO(3/68) -2 

80 

IZMIRAN(3/68) 

100 
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The GSFC model was given the 40y weight even though it had 
a 61y residual to the 0G0-4 data since the 0G0-2 figure was 39y, 
the other satellite residuals were low, and it has the overall 
lowest residuals to the surface data, The AFCRL model and RGO 
contributions were roughly equivalent but the AFCRL was given a 
slightly smaller figure due to its lower residual to 0G0-2, 

COSMOS -49 and the surface data. The IZMIRAN model was assigned 
a slightly higher weight because of its uncertainty in the polar 
regions due to its being derived from data at less than 50° latitude, 
This possible difficulty is evidenced by its relatively high 
residuals to data sets containing polar contributions (e.g., OGO-2, 
OGO-4, observatory, land/sea, Scandinavian airborne, and Project 
MAGNET). 

There was less basis for rational comparisons in combining 
the secular change terms. Hence each model previously used was 
weighted equally and the USCGS and RGO-1 models included since the 
secular change was independently derived for each. 

Although more lengthy considerations may have resulted in an 
improved procedure for deriving the first IGRF, this formulation 
provided a model composed of some contribution from each organization 
and at the same time, within the restrictions as to the number of 
coefficients, produced a model which agrees tolerably well with the 
test data set. This agreement is seen in the last column of Table 4. 
Surprisingly, the procedure appeared to produce a residual equal to 
or less than that of the contributing models for some of the data sets 


at the n* = 8 truncation level. 
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T he Resulting Model IGRF(lQ/68) 

Since the IGRF is a composite of several models, it can be 
compared with each as given in Table 5, Here is listed for each 
of the contributing coefficient sets the deviation from the 
resulting IGRF. Although the disagreements between the various 
terms are sometimes relatively large for those with amplitudes of 
the order of 1 to lOy, those of higher magnitude are surprisingly 
close. Of the main field terms there seems to be the largest 
discrepancy between those having m » 1, particularly for the 1ZMIRAN 
model with n « 1, 2, 3, 4, and 6. 

The final IGRF (10/68) coefficients are given in Table 6 and 
maps of the field and its secular change given in Appendix 1, 

Appendix 2 gives a possible minor modification based on a suggested 
change of scale to a standard mean earth radius of 6371 in place of 
6371.2 km. 

Recommendations 

We would like to make a few recommendations as to the way an 
international reference field might be used, As can be seen in this 
report and others we have published ( Cain et al., 1965 ; Cain et al. , 1967; 
Cain and Hendricks , 1968), ambient values of the earth's field are 
dependent on contributions from the core, crust, subsurface, and 
ionospheric electric currents, and the effects of trapped plasma, 
magnetospher ic boundary, and tail effects. The exact secular variation 
is subject to shifts which make a linear fit with time increasingly 
uncertain beyond a few years. Further, even for the decade of validity 
of the IGRF, 1960-1969, we already know that there are more accurate 


models available. 
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-3 

8 

8 

-3 

3 

X 

8 

6 

-9 

-16 


TABLE 5(b) 



OT 


n 

m 

GSFC ( l 2 / 66 ) 

AFCRL ( 3 / 68 ) 

RGO -2 

XZM 7 RAN 

US COS 

RGO -1 

1 GRF 

1 

0 

- 1,5 

- 3.3 

1,7 

3,3 

- 0.3 

0.2 

15,3 

1 

1 

0.7 

- 1.9 

0,3 

- 1,0 

1,8 

- 0,4 

8.7 

2 

0 

0.3 

«*3 , 3 

0,4 

- 0,7 

3,3 

- 2,2 

- 24.4 

2 

1 

- 0,5 

- 0.5 

- 1.3 

- 0,2 

2.4 

-1 , 6 

0,3 

2 

2 

- 4,2 

- 2.0 

3 .C 

0.4 

2,1 

2,9 

- 1,6 

3 

0 

- 1.7 

0.1 

- 0.2 

0,6 

1.1 

- 0.2 

0,2 

3 

1 

0,0 

0.6 

0.8 

1,9 

- 3.3 

1,3 

- 10,8 

3 

2 

1.7 

0,5 

- 1,7 

- 0.7 

0.2 

- 2.6 

0,7 

3 

3 

- 3.1 

- 3.0 

2.8 

2.1 

1.3 

3.2 

- 3.8 

4 

0 

2,1 

1.4 

- 0.3 

0,9 

- 4.1 

1,3 

- 0,7 

4 

1 

0.5 

0,6 

- 0.2 

0.1 

- 1,0 

0,8 

0,2 

4 

2 

1,3 

2.6 

- 1.0 

0.7 

- 3.6 

0.8 

- 3.0 

4 

3 

0.8 

- 0,4 

1.1 

0.1 

- 1.7 

0.3 

- 0.1 

4 

4 

- 1.4 

1.2 

- 0,9 

0.8 

0,3 

- 0,9 

- 2.1 

5 

0 

-0 , 0 

1,4 

- 0,9 

- 1.3 

0.7 

- 1.1 

1,9 

5 

1 

- 0,5 

- 0,5 

- 1.1 

- 0.8 

2.9 

- 0.7 

1.1 

5 

2 

0,9 

- 0,2 

0.1 

- 1.8 

0.9 

- 1.3 

2.9 

5 

3 

- 0,6 

- 0,2 

- 1,6 

- 1.0 

3.2 

- 0.9 

0.6 

5 

it 

- 0.6 

1.5 

- 1.0 

- 0,7 

0.8 

- 1.1 

0.0 

5 

5 

0.8 

2,7 

- 0.3 

- 1.2 

- 2,1 

0.4 

1.3 

6 

0 

- 0,4 

rl.l 

0,1 

- 0.4 

1.8 

0.1 

- 0.1 

6 

1 

1.2 

- 0.7 

0.3 

0.6 

- 1.3 

- 0,2 

- 0.3 

C 

2 

- 0.2 

- 0,6 

- 0.1 

- 0.4 

1.2 

0.7 

1.1 

6 

3 

0.7 

- 0.7 

0.1 

- 0.6 

0.4 

- 0,5 

1.9 

6 

4 

1.2 

0.9 

0.4 

0.8 

- 3.4 

1.0 

-0 , 4 

C 

5 

0.6 

- 2.3 

0.4 

0.4 

0.9 

0.6 

- 0.4 

6 

G 

0 . G 

0.6 

- 0.8 

0.2 

- 0.6 

- 2.0 

- 0,2 

7 

0 

- 0.1 

0.7 

0*5 

0.5 

- 1.6 

0.5 

- 0.5 

7 

X 

- 0.0 

- 0,0 

0.3 

0,3 

- 0,6 

0.3 

- 0.3 

7 

2 

- 0.9 

- 0.3 

0.7 

0,7 

- 0,2 

0.7 

- 0.7 

7 

3 

- 0.4 

1.3 

0.5 

0.5 

- 1.9 

0.5 

- 0.5 

7 

4 

- 0.3 

- 0.8 

- 0.3 

- 0.3 

1.5 

- 0.3 

0.3 

7 

5 

- 0.7 

- 0,7 

0.0 

0.0 

1.4 

0.0 

- 0.0 

7 

6 

0.0 

0.9 

0.2 

0.2 

- 1.3 

0.2 

- 0.2 

7 

7 

- 0.1 

- 1.6 

0.6 

0.6 

0.5 

0.6 

- 0.6 

8 

0 

0.3 

- 0.3 

- 0.1 

- 0.1 

0.3 

- 0,1 

0.1 

8 

1 

0.1 

- 0.1 

- 0,4 

- 0.4 

0.8 

- 0.4 

0.4 

8 

2 

1.3 

- 0.5 

- 0.6 

- 0.6 

0.3 

- 0,6 

0,6 

8 

3 

- 0.2 

- 0,8 

- 0.0 

- 0,0 

1.1 

- 0,0 

0.0 

8 

4 

0.5 

- 0,8 

0.0 

0.0 

0.2 

0.0 

- 0.0 

8 

5 

- 0.0 

- 0.1 

0.1 

0,1 

- 0.0 

0.1 

- 0.1 

8 

G 

0.2 

- 0.7 

- 0.3 

- 0,3 

1.1 

- 0.3 

0.3 

8 

7 

0.1 

- 1.4 

0.3 

0.3 

0.6 

0.3 

- 0.3 

8 

8 

0.0 

- 0,4 

0.5 

0,5 

- 0,6 

0.5 

- 0.5 


TABLE 5(c) 


HT 


n 

m 

GSFC ( 12 / 66 ) 

AFCRL < 3 / 68 ) 

1 

0 

0.0 

0.0 

1 

1 

- 1,7 

1.4 

2 

0 

0,0 

0.0 

2 

1 

- 2.3 

- 4.1 

2 

2 

0,0 

1.3 

3 

0 

0,0 

0.0 

3 

1 

1.5 

0.8 

3 

2 

1.8 

1.4 

3 

3 

1.1 

- 3.3 

4 

0 

0.0 

0.0 

4 

1 

- 2.1 

- 1.4 

4 

2 

- 1,4 

- 3.0 

4 

3 

- 1.2 

- 0.2 

4 

4 

- 2,6 

- 0.3 

5 

0 

0,0 

0.0 

5 

1 

- 0.3 

- 0.1 

5 

2 

- 0,1 

0.5 

5 

3 

- 0,3 

0,6 

5 

4 

- 0,3 

- 2.6 

5 

5 

0.1 

0.2 

6 

0 

0.0 

0.0 

6 

1 

1,1 

0.5 

6 

2 

0.4 

0.5 

6 

3 

0.6 

- 0,1 

6 

4 

- 0.3 

1.3 

6 

5 

Q , 4 

0,5 

6 

6 

- 0.1 

- 1.3 

7 

0 

0.0 

0.0 

7 

1 

0.1 

0.5 

7 

2 

- 0.2 

- 1.1 

7 

3 

0.1 

0.5 

7 

4 

0.6 

- 0 .G 

7 

5 

- 0.8 

- 0.8 

7 

G 

0.4 

1.5 

7 

7 

0.7 

1.1 

8 

0 

0,0 

0.0 

8 

1 

- 0.7 

- 0.4 

8 

2 

- 0.1 

- 0.0 

8 

3 

0.4 

- 0.7 

8 

4 

- 0.7 

- 0.2 

8 

5 

0.4 

- O.i 

8 

6 

0.5 

0,2 

8 

7 

- 0.1 

- 1.0 

8 

8 

- 0.2 

- 0.4 


RGO -2 

IZMIRAN 

USCGS 

RGO-J 

0.0 

0.0 

0.0 

0.0 

2,3 

- 0,2 

- 1.7 

2.9 

0,0 

0*0 

0.0 

0,0 

- 0.2 

2.0 

4.7 

0,4 

- 1 .3 

0.9 

- 1.1 

- 1.5 

0.0 

0.0 

0.0 

0.0 

0.8 

0.1 

- 3.1 

-1 . 0 

0.3 

0.0 

- 3.7 

0.9 

- 1.3 

2,1 

1.5 

- 0.8 

0.0 

0.0 

0.0 

0.0 

1.1 

1.5 

1.0 

3.1 

- 0.6 

- 0.7 

5.7 

- 2,3 

0,1 

- 0,4 

1.7 

- 0.2 

1.2 

2.6 

- 0.8 

1 • 5 

0.0 

0.0 

0.0 

0,0 

- 0.3 

- 1.9 

2.6 

- 0.4 

0.3 

-0 , 3 

- 0,3 

0,6 

0.4 

1.2 

- 1.9 

0.6 

0.2 

- 0.6 

3.4 

0.6 

0.3 

0.5 

- 1.1 

0.8 

0.0 

0.0 

0,0 

0.0 

- 1.1 

0.5 

- 1.0 

- 1.3 

1.4 

0.6 

- 2,9 

0 . 6 

- 0.0 

- 1.1 

0.8 

- 1.1 

0.1 

1.1 

- 2.1 

- 1.0 

- 0.1 

- 0.1 

- 0.7 

- 0,2 

0,1 

- 0.9 

2.3 

- 1.1 

0.0 

0,0 

0.0 

0,0 

1.1 

1.1 

- 2.7 

1.1 

- 0.3 

- 0.3 

1.9 

- 0,3 

- 0.4 

- 0.4 

0.2 

- 0,4 

- 0.2 

- 0.2 

0.4 

- 0.2 

- 0.4 

- 0.4 

2.4 

- 0.4 

- 0.2 

- 0,2 

- 1.4 

- 0.2 

- 0.3 

- 0.3 

- 1.2 

- 0,3 

0.0 

0.0 

0.0 

0.0 

- 0.1 

- 0.1 

1.3 

- 0,2 

0,2 

0.2 

- 0,2 

0.2 

0.3 

0.3 

- 0.4 

0.3 

0.2 

0.2 

0.5 

0.2 

0,3 

0.3 

- 1.0 

0.3 

0.4 

0.4 

- 1.5 

0.4 

0.3 

0.3 

0.5 

0.3 

0.3 

0.3 

0.1 

0.3 


IGRF 

0,0 

- 2.3 

0.0 

- 11.8 

- 16,7 

0.0 

4.2 
0.7 

- 7,7 

0.0 

- 0.1 

1.6 

2.9 

- 4.2 

0.0 

2.3 
1.7 

- 2.4 

0.8 

- 0.3 

0.0 

- 0.9 

- 0,4 

2.0 

- 1.1 

0.1 

0.9 

0.0 

- 1.1 

0.3 

0.4 

0.2 

0.4 

0,2 

0.3 

0.0 

0.1 

- 0.2 

- 0.3 

- 0.2 

- 0.3 

- 0.4 

- 0,3 

- 0.3 


TABLE 5(d) 



EPOCH - 1965.0 


I * G , R , F ,(10/68) 


N 

M 

G 

H 

GT 

HT 

1 

0 

-30339 

0 

15,3 

0.0 

1 

1 

-2123 

5758 

8.7 

- 2.3 

2 

0 

-1654 

0 

- 24.4 

0.0 

2 

1 

2994 

-2006 

0.3 

- 11.8 

2 

2 

1567 

130 

- 1.6 

- 16.7 

3 

0 

1297 

0 

0.2 

0.0 

3 

1 

-2036 

-403 

- 10.8 

4.2 

3 

2 

1289 

242 

0.7 

0.7 

3 

3 

843 

-176 

- 3.8 

- 7.7 

4 

0 

958 

0 

- 0.7 

0.0 

4 

1 

805 

149 

0.2 

- 0.1 

4 

2 

492 

-280 

- 3.0 

1.6 

4 

3 

-392 

8 

- 0.1 

2.9 

4 

4 

256 

-265 

- 2.1 

- 4.2 

5 

0 

-223 

0 

1.9 

0.0 

5 

1 

35 7 

16 

1.1 

2.3 

5 

2 

246 

125 

2.9 

1.7 

5 

3 

-26 

-123 

0.6 

- 2.4 

5 

4 

-161 

-107 

0.0 

0.8 

5 

5 

-51 

77 

1.3 

- 0.3 

6 

0 

47 

0 

- 0.1 

0.0 

6 

1 

60 

-14 

- 0.3 

- 0.9 

6 

2 

4 

106 

1.1 

- 0.4 

6 

3 

-229 

68 

1.9 

2.0 

6 

4 

3 

-32 

- 0.4 

- 1.1 

6 

5 

-4 

-10 

- 0.4 

0.1 

6 

6 

-112 

-13 

- 0,2 

0.9 

7 

0 

71 

0 

- 0.5 

0,0 

7 

1 

-54 

-57 

- 0.3 

- 1.1 

7 

2 

0 

-27 

- 0.7 

0.3 

7 

3 

12 

-8 

- 0.5 

0.4 

7 

4 

-25 

9 

0 . 3 

0.2 

7 

5 

-9 

23 

- 0.0 

0.4 

7 

6 

13 

-19 

- 0.2 

0.2 

7 

7 

-2 

-17 

- 0.6 

0.3 

8 

0 

10 

0 

0.1 

0.0 

8 

1 

9 

3 

0.4 

0.1 

8 

2 

-3 

-13 

0.6 

- 0.2 

8 

3 

-12 

5 

0.0 

- 0.3 

8 

4 

-4 

-17 

- 0.0 

- 0.2 

8 

5 

7 

4 

- 0.1 

- 0.3 

8 

6 

-5 

22 

0.3 

- 0.4 

8 

7 

12 

-3 

- 0.3 

- 0.3 

8 

8 

6 

-16 

- 0,5 

- 0.3 


TABLE 6 
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The IGRF was developed as a result of the request of those who 
would wish a standard field model where the permanence of a standard 
over a period of years outweighs the advantages of a high accuracy. 
Thus the ultimate use of this model and further requests for revisions 
must be left to the users. 

The way to test whether IGRF (10/68) is suitable to any particular 
need is to periodically test newer or more accurate models and to 
decide on the basis of the differences whether the continued use is 
adequate. As the core field deviates more and more from the IGRF 
estimate, the accuracy will continuously decrease. 

We have already made this test in regard to analysis of the time 
variations of the COSMOS -49, 0G0~2, and OGO-4 data. For such studies 
the IGRF is quite useless, even the GSFC (12/66) model is insufficient, 

l 

and fits based on the data themselves are being used. For higher 
accuracy studies, we would suggest using the GSFC(12/66) model over 
the range 1900-1965 and the POGO (10/68) model from 1965/ through 1968. 
Beyond 1968, POGO(10/68) could be used until it is updated by more 
recent data and planned improvements in the analysis. 

The computations of the magnetic field from the IGRF or other 
magnetic field coefficients can be effected using a wide variety of 
computer programs currently available. One such set of programs 
based on a code originally developed by Jensen and Whitaker (1960), 
may be obtained from: 

World Data Center A for Rockets and Satellites 
Goddard Space Flight Center (601) 

Greenbelt, Maryland 20771 
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These codes convert the Schmidt normalized coefficients 
internally in the computer to a more efficient Gauss normalized 
form, update them to the epoch requested, and compute the 
geocentric components from the scalar gradient of the potential 
function given the geocentric position. Conversions are also 
provided so that one can enter the programs with a geodetic 
position which is then converted to geocentric, and also for 
rotating the output geocentric components into geodetic directions. 
Ignoring the differences between geodetic and geocentric coordinates 
will create errors up to about 200y. 



APPENDIX I 


Main Field Component and Isoporic Charts 
Computed From IGRF (10/68) for 1965.0 at the Earth's Surface 


The following figures represent the surface contours of the 
various geodetic components of the geomagnetic field and its 
secular change as computed by the IGRF . These diagrams are very 
similar to those given by Cain and Hendricks (1968) for the 
GSFC(12/66) field and are drawn automatically using a computer 
program originally used for weather maps ( Cain and Neilon , 1963). 

The plots thus are drawn including the algebraic "lows” and 
"highs” of the component being displayed. These extrema occur at 
the center of the or symbols. The dip poles are noted for 


the H chart as 
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IGRF 1965.0 
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IGRF 1965 
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IGRF 1965 












SECULAR CHANGE OF TOTAL FORCE (gammas/year) 



IGRF 1965 









APPENDIX 2 


Coefficient Normalization 


All of the previous field derivations have arbitrarily set 
the earth’s mean radius at 6371.2 for the value of a in the 
factors (a/r) n+ '*' of the potential expansion. This value stemmed 
from the old standard earth constants with equatorial radius 
6378.388 and flattening 1/297. However, the new constants have 
become 6378.15 and 1/298.25 respectively. Integrating 


we obtain 


F 

F 



cos© d0 


~ ln(m + a/b) 


where 


m = JaP - b s /b 

a = equatorial radius, and 

b = a(l - f) is the polar radius 
with f the flattening factor 


The values with the old and new constants are as follows: 


f a b F 

297 6378.39 6356.91 6371.21 

298.25 


6378.16 


6356.77 


6371.02 




APPENDIX 2 
(cont ’d) 


I would like to recommend that for the sake of simplicity 
and to assure that we are not bound to constants of only historical 
significance, that we adopt the value of 6371 for a. This is a 
very slight change and has the effect of only altering the gj term 
from -30339 to -30342 and the hj term from 5758 to 5759. The 
constants a to make the correction g = g' + ag' where g' is the 
old values of g or h, are 

n <*X10 B 

1 9 


2 

13 

3 

?6 

4 

19 

5 

22 

6 

25 

7 

28 

8 

31 

9 

35 

10 

38 

11 

41 

12 

44 
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